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Abstract 

This article concerns the correspondence between thermodynamics and the mor- 
phology of simple fluids in terms of clusters. Definitions of clusters providing a 
geometric interpretation of the liquid-gas phase transition are reviewed with an eye 
to establishing their physical relevance. The author emphasizes their main features 
and basic hypotheses, and shows how these definitions lead to a recent approach 
based on self-bound clusters. Although theoretical, this tutorial review is also ad- 
dressed to readers interested in experimental aspects of clustering in simple fluids. 
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1 Introduction 



Simple fluids are classical systems in which chemically inert particles interact 
through a pairwise potential. The main features of the potential are the hard- 
core repulsion and the short-ranged attraction which give rise to thermal phase 
transitions. For example, noble gases and alkali metals are often modeled as 
simple fluids. 

Moreover, attraction between particles promotes cluster formation. It is then 
sensible to look for a relation between thermodynamics and the morphology 
of a fluid in terms of clusters. In particular, a geometric interpretation of 
the liquid-gas phase transition is an important issue in statistical mechanics. 
To put it more precisely, there are two main issues: firstly to understand 
condensation as the sudden formation of a macroscopic cluster, and secondly 
to describe the morphology of the fluid at the critical point with a view to 
infer its thermodynamic properties. 

To begin, it must be said that a geometric description of a fluid is outside the 
scope of thermodynamics. In addition to standard tools of thermal statistical 
physics, one has to define what is a cluster. At first sight, it may seem obvious 
to define a cluster as a set of particles close to each other and far from other 
clusters. However, this definition becomes quite ambiguous at high density, 
even below the critical density. Therefore, more sophisticated definitions of 
clusters were introduced, giving the geometric description of simple fluids a 
long and rich story. 

This article presents a comprehensive review of proposed definitions of clusters, 
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with a view to establish a relation between thermodynamics and a geometric 
description of the fluid. However, the reader will see that this purpose will 
take us beyond this particular topic. On the one hand, our aim is to out- 
line the main features of these various clusters. In particular, we highlight 
the hypotheses on which their definitions are built, in order to discuss their 
physical relevance. Indeed, some of these definitions were designed as mathe- 
matical tools to study thermodynamic properties. However, these clusters are 
often used in literature as though one could observe them experimentally. On 
the other hand, we clarify the relations between these definitions by follow- 
ing a chronological order so as to avoid frequent misunderstandings. Hence, 
we develop a new point of view about the Kertesz line associated with the 
Coniglio-Klein clusters. From this standpoint, we present a recent approach 
in terms of self-bound clusters that provides a physical interpretation of the 
liquid-gas phase transition and suggests the existence of a percolation line in 
the supercritical phase of simple fluids. 

This article is aimed not only at readers interested in fundamental aspects of 
the relation between thermodynamics and geometry, but also at experimental- 
ists who deal with clusters in various fields of physics. In particular, we think 
of physicists who search for vestiges of phase transitions in small systems like 
atomic nuclei, and aggregates. 

The plan of the review is as follows. In section 2, we recall the definition of 
a simple fluid and describe its typical phase diagram. Section 3 is devoted to 
definitions of clusters based on the "perfect gas of clusters model" developed 
in appendix A. Section 4 discusses microscopic definitions of clusters proposed 
in the framework of the lattice-gas model by means of percolation theory. In 
section 5, we are interested in self-bound clusters defined by energetic criteria. 
We present the correspondence between the phase diagram and their cluster 
size distribution, in the framework of the lattice-gas model, and by using 
numerical simulations of a Lennard- Jones fluid. Finally, section 6 contains a 
summary and we close by discussing several open questions. 



2 From thermodynamics to the morphology of simple fluids 



A simple fluidQ is a classical system composed of N structureless particles 
of mass m, interacting through a two-body additive central potential u(rij), 
where is the distance between particles i and j. The potential must have 
a repulsive hard-core and a short range attraction that vanishes faster than 
—l/rfj to ensure the thermodynamic limit to exist in three dimensions ( |Kucllcj , 



For a comprehensive review about simple liquids, we refer the reader to the book 



of Hansen and McDonald ( 1986 ). 
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1969Q . The Lennard- Jones potential, for instance, satisfies these conditions: 
«r«) = M(^) 12 - (f ) 6 ] (1) 

where the two constants e and a fix the energy and length scales respectively. 
The Hamiltonian of the fluid is the sum of a kinetic and an interaction energy 
term: 

N v 2 

+ £«fa) (2) 
i=i zm i<j 



where fi and pi are the position and momentum coordinates of the i particle, 
and Tij = \fi — fj\. In the canonical ensemble, for a given volume V and 
temperature T, the partition function of the fluid is given by 

Qn ( T > V) = iVlW / e " W,},{ft}) ^ • • • d ^dpi ■ ■ ■ dp N 

where f3 = 1/kT, k is the Boltzmann's constant and h the Planck's constant. 
Integrations over the momentum coordinates are straightforward and allow us 
to write the partition function as 

Qn(T,V) = -± w Z n (T,V) (3) 



where 



\j2iirakT ^ 



and 



Z N (T, V) = [ IJ e-^dn . . . dr N . (5) 

V {< 3 



Equation (4) defines the thermal wavelength. The thermodynamics of the fluid 
is then determined by the behaviour of the configuration integral Z^(T, V) in 
the thermodynamic limit. 

The typical phase diagram of the system is shown in Fig. 1 . In the P — T dia- 
gram, coexistence curves separate the plane into three regions corresponding 
to the solid, liquid, and gas states. The liquid-gas coexistence curve stops at 
the critical point (C). Crossing this curve, the system undergoes a first order 
phase transition signaled by a latent heat and the coexistence of a low density 
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phase, the gas, and a high density phase, the liquid. In other words, compress- 
ing a gas at constant temperature T below the critical temperature T c , a first 
drop of liquid appears when, in the phase diagram, one reaches the liquid- 
gas coexistence curve at the pressure P con d(T), at the corresponding density 
Pcond{T). This phenomenon is called "condensation". On the other hand, just 
at the critical point, a continuous phase transition occurs, while the compress- 
ibility and the specific heat diverge. The critical behaviour is known to be in 
the universality class of the Ising model (see for instance |GoldenfeId| , |1992| ). In 
the T — p diagram, above the bell-shaped liquid-gas coexistence curve, there 
is one single fluid phase called the "supercritical phase" .0 Because the ther- 
modynamic properties of the system vary smoothly along any path that does 
not cross a coexistence curve, it is possible to pass continuously from the gas 
to the liquid phase by following a path such as the dashed line shown in Fig. 
1. 
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Fig. 1. Schematic phase diagram of a simple fluid in the P — T plan (left) and the 
T — p plan (right), where P, T and p are respectively the pressure, temperature and 
density of the system. The solid, liquid and gas phases are respectively noted s, I, 
and g. The critical point is C. The dashed line with the arrow represents a path in 
the phase diagram that does not lead to singularities of the partition function. 

In the gas, liquid, and supercritical phases, the system is homogeneous, density 
is uniform on average. However, attractive interaction between particles con- 
tributes to forming localized morphological structures, the clusters. It must be 
emphasized that in order to study the morphology of a fluid, we need to give 
a definition of clusters in addition to standard thermodynamic tools. Indeed, 
the partition function depends only on the sum of the two-particle interac- 
tions, without discerning if a particle belongs to a cluster or not. In other 
words, one could arbitrary distribute the N particles among a set of clusters, 
while the partition function and all the thermodynamic properties would not 
be altered. Of course, the choice of the definition of clusters has no influence 
on the thermodynamic behaviour of the fluid, but depends on the problem we 
want to study. 



The supercritical phase of a fluid (water) was discovered by the Baron Cagniard 



de La Tour (|1822|) . 
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Historically, the first motivation for studying the morphology of a fluid was to 
interpret condensation as the sudden formation of a macroscopic cluster. At 
that time, the purpose was especially to determine thermodynamic quantities 
with the help of "effective" clusters used as a convenient and powerful tool. 
These first definitions of clusters are based on the "perfect gas of clusters 
model" . 



3 Perfect gas of clusters 



The basic idea is to assume that an imperfect fluid made up of interacting 
particles can be considered as a perfect gas of clusters at thermodynamic and 
chemical equilibrium: clusters do not interact with each other and do not 
have any volume. We shall discuss the validity of the perfect gas of clusters 
approximation at the microscopic level when we tackle the problem of self- 
bound clusters in section 5.2. 

Let us suppose we deal with a collection of non-interacting clusters, without 
knowing how we can build them as sets of interacting particles. Clusters of 
given size s are characterized by their mass ms, a chemical potential fi s , and 
a partition function q s (T, V). We show in appendix A that the pressure, the 
density, and the cluster size distribution, that is the mean number of clusters 
of size s, are respectively given by 



1 oo 

P p =^E9/ (6) 

V s=l 

oo 

P =yT, S( ls^ (7) 
s=l 

n s (T,z) = q s z s (8) 

where z = e 13 ^ is the fugacity and /i the particle chemical potential. The 
equation of state of the fluid can be obtained by eliminating z in the pressure 
expression with the help of Eq. (7). It is interesting to write Eq. (6) as a 
function of the cluster size distribution given by Eq. (8): 

oo 

(3PV = Y / n s . (9) 



Equation (9) shows strikingly that the imperfect fluid can be seen as a perfect 
gas composed of m = J2^Li n s non-interacting clusters. 

According to Eqs. (6)- (8), thermodynamic quantities, as well as the cluster 
size distribution, are completely determined by the partition function of the 
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clusters. Now, two points of view can be considered. In the pioneer theory of 
condensation proposed by Mayer and his collaborators in 1937, clusters result 
from an exact enumeration of all the possible combinations of interacting 
particles (sec. 3.1). It is equivalent to giving a formal expression of q s (V,T) 
that is impracticable to handle for large size clusters. On the other hand, 
Frenkel, Band (sec. 3.2), and Fisher (sec. 3.3) propose a phenomenological 
expression of the partition function q s (V,T), and assume that clusters do not 
interact whatever the density. 



3.1 Mayer clusters 



Mayer's theory of condensation is described in many textbooks of statistical 
physics (see for example |Pathria , 1972 ; Huang , 1987 ) and in Mayer's arti- 
cles ( |Mayer| , |1937| ; [Mayer and Ackermannj , |1937| ; [Mayer and Harrison| , |1938| ; 
Harrison and Mayeij , |1938| ; Mayer and Mayer] , |1940|) . Here, we just want to 
show how Mayer clusters can be seen as a particular case of the perfect gas of 
clusters model. 



Originally, this theory was intended to express the exact equation of state of a 
real fluid, made up of interacting particles, as a series expansion in the density. 
To this end, Mayer decomposes the partition function (3) into the sum over 
all the possible partitions of particles into independent mathematical clusters. 
The first step is to introduce a function of the potential that is significant only 
if particles are close to each other: 

fa = e-Mr*) _ i 

The function \fa\ is bounded everywhere, unlike the potential u(ry), and 
tends to when becomes large compared to the interaction range. Using 
this function as a small parameter, one can write the configuration integral 
(5) as a diagrammatic expansion by associating an iV-particle graph with each 
term.0 A graph consists of vertices and bonds: to each particle corresponds a 
vertex, and the function is represented by a bond connecting vertices % and 
j. A set of pair by pair connecting vertices is called a "Mayer cluster". Each 
graph, or term of the expansion, is characterized by a cluster size distribution. 

By collecting the terms having the same cluster size distribution, one obtain a 
new expression of the partition function that can be simplified by moving over 
to the grand canonical ensemble. In the thermodynamic limit, pressure and 
density are eventually written as series expansion in the fugacity. As a result, 
calculations lead to the expressions (6) and (7), in which q s /V is replaced by 

4 According to Huang ( |1987j ), this is the first graphical representation of a pertur- 
bation series expansion in physics. 
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b s (T)/X 3 , where coefficients b s (T) are the so-called "cluster integrals" intro- 
duced by Mayer, and A is the thermal wavelength of the particle defined by 
Eq. (4). 

In order to locate the condensation point, that is the density at which a 
macroscopic cluster appears, we have to examine the cluster size distribution 
given by Eq. (8). At a given temperature, the behaviour of the coefficients 
b s (T) for s > > 1 shows that with increasing density, n s (T, p) becomes non- 
zero when p > pcond- A macroscopic cluster appears, and Mayer identifies this 
particular density p con d with the density of condensation. However, Yang and 
Lee ( |1952| ) proved that Mayer's theory is exact in the gas phase, for p < p CO nd, 
but cannot be carried forward into the liquid phase. 



By definition, Mayer clusters form a perfect gas of clusters and are designed 
to evaluate thermodynamic quantities. Yet, it can be shown that b s (T), and 
therefore the cluster size distribution, may be negative in certain thermody- 
namic conditions. There then appears to be no physical interpretation of these 
clusters (FTsEerl |1967a|; |Pathria|, |1972|). 



3.2 Frenkel clusters 



In contrast to Mayer's theory, the basic idea of Frenkel's modelQ ( ]F renkel , 
I19394 R) is to use directly the perfect gas of cluster model by choosing a 
phenomenological expression of the partition function of clusters. Afterward, 
this model was modified to allow for the volume of the clusters (|Band| , |1939b| ; 



Stillingci] , |1963|) , but attraction between clusters was still overlooked. 



By disregarding the degrees of freedom associated with the particles, Frenkel 
assumes the clusters of size s to be compact and writes their potential energy 
Ep s , for s sufficiently large, as the sum of a bulk and a surface term 



Ep s (T) = -e v s + e a s 



for 



s » 1 



where e v > and e a > are the bulk and surface potential energy by particle 
respectively. Frenkel does not take into account the entropy of the clusters and 
infers the partition function by integrating over the position and momentum 
coordinates of the center of mass of the cluster: 



q s (T,V) = V( 



A ' 



-liEp a 



A^ 



5 A similar model has been independently proposed by Bijl (1935) and Band 
(|1939a|) . 
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where A/a/s is the thermal wavelength of a cluster of size s (see Eq. (4)). 
According to Eq. (8), the Frenkel cluster size distribution is given by 



n s {T,z) s3 



3 

2 



v T^ x < 10 > 

with 

y(T, z) = ze^ and x(T) = e~ Pea . 

The function y depends on the temperature and the density through z. The 
function x is independent of the density and always less than 1. For a given 
temperature T, x is fixed, the cluster size distribution depends on density only 
through y(T,z). When y < 1, n s decreases exponentially for s » 1: there 
is no macroscopic cluster, this corresponds to the gas phase. For y > 1, n s 
decreases as long as s is smaller than a particular size, and then increases 
exponentially: a macroscopic cluster appears, which indicates the formation 
of the liquid phase. The density at the condensation point is then p = p con d(T) 
such that y = y cond = 1. 

The results of Frenkel's phenomenological model are for the most part the 
same as those of the Mayer's theory of condensation. The main advantage 
is that Frenkel clusters do not have the pathology of the Mayer clusters: the 
cluster size distribution (10) is positive whatever the temperature and density. 
However, this model does not allow one to locate the critical point and describe 
the morphology of the fluid in this state. 

3.3 Fisher droplets 



As Frenkel's model, Fisher's droplet model ( [Fishcr| , |1967a| Jo|, |1971|) is based on 



the perfect gas of clusters approximation. However, Fisher writes the partition 
function of the clusters of size s with the help of additional features: he allows 
for the entropy of the clusters, clusters are not assumed to be compact, and a 
corrective term varying like Ins is added to the free energy of the clusters of 
size s, in addition to the surface and volume terms. 



Here, we shall not follow Fisher's original calculations (|Fisher|, |1967a|). We 



rather propose a simplified and more direct way of calculating the partition 
function of clusters of size s (|Fisher| , |1971| ). 

Let us first write the free energy of the clusters of size s. The mean internal 
energy and entropy of a cluster of size s >> 1, with a mean surface area A s , 
are written as the sum of a surface and volume term: 



U s = -u v s + u a A s 
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S s S V S ~\~ S a j\- S 

where u v > and s v > are the volume energy and entropy by particle 
respectively, and u a > and s a > 0, the surface energy and entropy by 
particle. Fisher does not assume the clusters to be compact and introduces a 
parameter a to characterize the mean surface of the clusters of size s: 

A S (T) = a {T)s a with < o < 1. 

The surface of a three-dimensional cluster should lie between the surface of a 
compact object (er = 2/3) and the surface of a chain (a = 1). Values of a less 
than 2/3 are interpreted by Fisher as an effective evaluation of the interaction 
between clusters. 

Fisher also adds to the free energy a corrective logarithmic term which has 
no obvious physical interpretation^ ( |Binder| , |1976b|) . The weight of this term 



is given by a second parameter r. We shall see later that this term plays a 
crucial role in this model. The free energy F S (T, V) of a cluster of size s is 
then given by 

—j3F s (T, V) = (3(u v + s v T)s — j3a (u a — s a T)s a — r Ins + lnc K 

The term proportional to In V results from the integration over the position of 
the center of mass of the cluster and cq is a constant. The partition function 
of a cluster of size s is eventually given by 

q a (T, V) = e -^(T,V) = c ^y s -t e f3{u v +s v T)s-f3a Q (u a -s a T)s° _ 

On the basis of the "perfect gas of clusters model" presented in appendix A, 
the pressure, density, and the cluster size distribution are written as 



1 oo oo 

PP = tt(T, z) = - £ q s z s = c V s x s °s~ T (11) 

V 8=1 S=l 

1 oo oo 

p=-Y / sq s z s = co^sy'x^s^ (12) 

V s=l s=l 

y = Y z =c y x s (13) 



where 



v(T,z) 
x(T) 



■ ze 



P(u v +s v T) 
-/3a (u a -s a T) 



6 This term is supposed to allow for the surface undulation of the clusters flEssam 
and Fisheij |1963| ; |Fisherj |1971|) . 
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Following Fisher, we have introduced the series ir(T, z) in Eq. (11) and we 
define the derivative of order n by 

a oo 

7r^(T, z) = (z— )<»>7r(T,z) = c J2y s x s °s n -\ 

The condensation point is determined by the appearance of a macroscopic 
cluster. According to the value of x, two cases are possible: if x < 1, as 
in Frenkel's model, the density of the condensation point p con d is given by 
y = Vcond = 1, that is z con d = e -f 3 ( u ^+^ T \ On the other hand, if x > 1, 
the cluster size distribution increases when y > 1 and the series (11) and 
(12), which give the pressure and density, diverge. Therefore, condensation 
only happens when x < 1, that is for T < T c = u a /s a . This upper limit 
on the condensation temperature is interpreted as the critical temperature. 
Consequently, Fisher's model is not valid for T > T c or p > p con d (^ > 1 or 

y>i)- 

Despite the basic approximation of the perfect gas of clusters model, Fisher 
extends the validity of his model to the case of high densities and studies the 
vicinity of the critical point. At the critical point, y — x — 1, the critical 
density is given by 

oo 

p c = 7T (1) (T C , Z cond ) = C sl ~ T - 

s=l 

The convergence of this series implies r > 2. The behaviour of the thermo- 
dynamic quantities in the neighbourhood of the critical point is obtained by 
moving over to the continuous limit. The series 7r^(T, z) can be written for 

% %cond- 

oo °9 

K {n \T,z cond ) = c^s n - T x sa ~ c / s n ~ T e- es °ds 

s=i i 

where 9 = lna; -1 = f3aos a {T c — T). Performing the change of variable t = #s CT , 
we have 

7r<">(r, w) - ^^r(^±l) ~ (T C -T)^P\ 

a o 

We show in appendix A how thermodynamic quantities are linked to the 
derivatives ir^ . According to the value of n, we infer the behaviour in the 
vicinity of the critical point of the specific heat C v (n — 0), the density (n— 1), 
the compressibility x ( n — 2), and the pressure (n = 0): 
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C v ~ (T c 




with 


«T = 


o 

z 


r 


Pc 


- P ~ (T c 


_ J^Ar 


with 


#r = 


r 


- 2 




(7 








with 


7t = 


3 


— r 






a 


Pc- 


-P~(p c 


-p) 5t 


with 


<5t = 




1 


r 


- 2 



- 1 



Fisher's droplet model allows one to write the thermodynamic critical expo- 
nents Qy, /3t, 7t, and 5^ as functions of the two geometric parameters a and 
r. Furthermore, the scaling relations between the thermal exponents ( |Fisherj , 
1967b| ), which had made this model well-known, follow easily from Eqs. (14): 



a T + 2/5 T + 7T = 2 (15) 
P T (8 T - 1) = 7T- (16) 

Let us now show the relation between the thermodynamics and the geometric 
description of the fluid at the critical point. From Eq. (13), the cluster size 
distribution is written for x = y = 1 

y = c s . 

We now see that the term in In s, added by Fisher to the free energy, is essential 
at the critical point. Without this term, we cannot infer the scaling relations 
(15)-(16) and the power law behaviour of the cluster size distribution. It has 
to be pointed out that the relation between a power law cluster size distribu- 
tion and the critical opalescence, observed at the critical point, is not clear 
today. A cautionary remark: this phenomenon, due to the divergence of the 
density fluctuations, is explained by thermodynamics without any reference 
to a geometric description of the fluid ( [Pathria| , |1972| ; IChandler], |1987j ). 



The two geometric exponents r and a can be written as functions of the 
thermodynamic critical exponents by reversing Eqs. (14): 



1 1 

a ■ 



/3t$t Pt + 7t 



St Pt + It 

In order to respect the critical behaviour of the thermodynamic quantities, 
exponents defined by Eqs. (14) must be positive. Consequently, the values of 
the geometric exponents are restricted to 2 < r < 3, and a > 0.5. Moreover, 
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the inflection pointQ of the critical isotherm implies 5t > 2. Therefore, 2 < 
r < 2.5, whatever the dimensionality ( [Kiang| , 1970Q . The exponents r and 



a, calculated from the values of the thermodynamic exponents, are given in 
Table 1. 
Table 1 

Values of the thermodynamic critical exponents and of the Fisher exponents a and 
t at the critical point, in mean field theory and in the universality classes of the 







0T 


IT 


St 


a 


T 


Mean field 





1/2 


1 


3 


2/3 


7/3 


2d— Ising 





1/8 


7/4 


15 


8/15 ~ 0.53 


31/15 ~ 2.07 


3d— Ising 


0.12 


0.325 


1.24 


4.8 


0.64 


2.21 



The link between thermodynamics and a geometric description of the fluid 
stems from the perfect gas of clusters model. Indeed, Eqs. (A.6)-(A.8) of ap- 
pendix A show that the first moments of the cluster size distribution are 
equal to thermodynamic quantities. For instance, moments of order 0, 1, and 
2 behave respectively as the pressure, the density, and the compressibility. 
Therefore, the second moment of the droplet size distribution diverges like 
the compressibility at the critical point. 

To test whether the droplet model is consistent with experimental data, the 
exponent r can be indirectly evaluated by means of the compressibility factor 
PP/p, measured at the critical point ( [Kiang| , |1970| ). Indeed, from Eqs. (11) 
and (12), we obtain for x = y = 1 



El 



(17) 



For a large variety of real fluids, the value of r inferred from Eq. (17) is in 
the order of 2.2, like the theoretical value calculated from the thermodynamic 
exponents (see Table 1). Moreover, a three-parameter fit of the compressibility 
factor along the coexistence curve provides a 1% agreement with experiment 
in the gas phase for water and carbon dioxide QRathjen et al.j |1972| ). 



However, we have to point out that Fisher's model has some important weak- 
nesses: 

• Fisher droplets are "effective clusters" which form a perfect gas of clusters. 



7 In the vicinity of the critical point, P c — P ~ (p c —p) ST , and the second derivative of 
the pressure as a function of the density must vanish, then d 2 P/dp 2 ~ {p c — p) St 2 = 
0. As a result, 5t > 2. 

8 However, Stauffer and Kiang ( |1971| ) extended the droplet model by taking into 
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This basic assumption casts some doubt upon the adequacy of the model 
at high density. Besides, as far as we know, no attempt as been made to 
observe the Fisher droplets in a real system. 

The equation of state produces some singularities in the supercritical phase, 
which is forbidden by the theorems of Yang and Lee ( |1952| ) ( [Kiang ancl 
Staufferl , pTCp-H 



The model does not describe the liquid phase beyond the coexistence curve 
for p > p con d ( [Kiang and Stauffer| , |1970| ; |Stauffer et al\ , |1971|) . 



Finally, despite these inadequacies, Fisher's phenomenological model provides 
a first morphological interpretation of the critical point in terms of clusters, as 
well as a thermodynamic study of the condensation curve based on geometric 
arguments. During the past three decades, Fisher droplets have become a guide 
for defining physical clusters. Indeed, we have to identify clusters among the 
particles of the fluid, in order to compare theoretical results with experiments 
and numerical simulations. In brief, we need a microscopic definition of clusters 
that behave like Fisher droplets at the critical point. 



4 Microscopic definitions of clusters in the lattice-gas model 



The power law behaviour of the Fisher droplet size distribution and the diver- 
gence of the mean droplet size at the critical point remind those of a percola- 
tion phase transition.^] Reader not familiar with percolation theory is referred 
to the book by Stauffer and Aharony (|1994 ). Hence, the works that followed 
those of Fisher have searched for a microscopic definition of clusters such that 
the thermodynamic critical point is a percolation threshold characterized by 
the Ising critical exponents. In other words, the microscopic clusters must 
verify the four following criteria at the critical point (pinder| , |1976a| ; [Kerteszj , 
gggg ; Etauft'erl , [OT]: |Staufl'er and Aharonyj gggg ): 



(1) The mean cluster size distribution behaves as a power law: n s ~ s~ T 
where r is the Fisher exponent (see Table 1). 

(2) A percolating cluster appears and its size S max varies like the order pa- 
rameter of the thermal phase transition, i.e the density difference between 
liquid and gas (or the spontaneous magnetization in the Ising model), 
with the exponent Pt- S max ~ (T c — T)^ T for T < T c . 



account the excluded-volume and an additional attraction between the clusters as 
a small perturbation. Critical exponents are not affected by these corrections. 

9 It is possible to cancel out these singularities by adding some corrective terms, 
without physical meaning, to the free energy of the clusters flReatto , 1970 ). 

10 Stillinger Jr ( 1963 ) put forward the idea for understanding the condensation as 
a percolation phase transition. In contrast, Fisher does not mention percolation in 
his original article ( Fisher , 1967a| ). 
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(3) The mean cluster size, i.e the second moment of the cluster size distribu- 
tion^ S mean ~ J2 S s2n si diverges as the compressibility (or susceptibility) 
with the same exponent 77-: S mean ~ \T — T C |~ 7T . 

(4) The connectedness length diverges as the thermal correlation length with 
the same exponent vt- 



We recall that in random percolation theory, the size of the largest cluster, 
the mean cluster size, the connectedness length, and the cluster size distribu- 
tion are respectively described at the percolation threshold by the following 
critical exponents /3 P , j p , u p , a p and r p . These exponents are universal for a 
given dimensionality and do not depend on the particular process of percola- 
tion (site/bond percolation, lattice/continuum space) flStaufier and Aharony) , 
1994| ). For the sake of comparison, the thermodynamic and random percola- 
tion exponents are respectively given in Tables 1 and 2. 



Table 2 

Values of the critical exponents in mean field theory and in the universality classes 
of random percolation theory in 2 and 3 dimensions (from Stauffer and Aharony\ , 



1994). 





Pp 


7 P 


O p 


Tp 


Mean field 


1 


1 


1/2 


5/2 


2d— Percolation 


0.139 


2.39 


0.39 


2.05 


3d— Percolation 


0.41 


1.80 


0.45 


2.18 



The simplest and most natural way for defining and identifying clusters among 
a group of particles is to base the criterion on proximity in configurational 
space ( |Band| , |1939a| Fisherj , |1971| ): two particles are linked if the distance 
between them is shorter than a given cutoff distance b. A cluster is then 
a group of particles linked pair by paiif* 2 ] (see Fig. 2). In this article, such 
clusters will be denoted "configurational clusters" . We now have to choose the 
value of the cutoff distance. At very low density, the identification of clusters 
depends little on it. In contrast, this definition becomes very sensitive to the 
choice of b at higher densities ( |Stillingerj , |1963| ; |Lee et aT\ , \1973[ |Abraham and 



Barker| , |1975j ; pinder| , |1975| ). As an illustration, if we were to double the cutoff 
distance in Fig. 2, all the particles would belong to the same cluster. 



11 In percolation theory, the mean cluster size is usually defined as milm\ where 
?7ifc = s k n s is the moment of order k of the cluster size distribution n s without 
the contribution of the largest cluster ( [Stauffer and Aharony , 1994). 

12 In the same spirit, a cluster can be defined either as particles lying inside a 
spherical shell centered on the center of mass of the cluster ( Reiss et al , 1968 ; Lee 
et al., 1973; Abraham, 1974; Abraham and Barker, 1975| ), or as a density fluctuation 



exceeding a given threshold (|Reiss et al\ , |1990[) . 



15 




Fig. 2. Particles in a box, in two dimensions. The radius of the discs centered at the 
positions of the particles is equal to b/2, where b is the cutoff distance. Configura- 
tional clusters are defined by the overlap of these discs. Hence, in this figure, one 
can see two monomers, three dimers and a cluster of size five. 



In order to apply the perfect gas of clusters model to a simple fluid, we could 
choose b ~ constant x r min , where r min is the distance corresponding to the 
minimum of the potential (|Binder| , |1975|) . Nevertheless, when a realistic po- 
tential is used, the attractive interaction between particles is never equal to 
zero, and there is no natural cutoff distance to defining the clusters: the choice 
of b is arbitrary. 



Thanks to its simplicity, the lattice-gas (Ising) model is a suitable framework 
to study configurational clusters and to test the four criteria of the Fisher 
droplets at the critical point. As we shall see in the next section, by choosing 
an attraction only between particles occupying nearest neighbour sites of a 
lattice, it becomes natural to choose a cutoff distance equal to the finite range 
of the interaction ( [Stoll et al.\ , |1972| ; |Binder| , |1975| ; [Miiller-Krumbhaar and 
~StoH , p76|) . 



4-1 Ising clusters 



First, let us recall the definition of the lattice-gas model. The M sites of 
a lattice are either empty or occupied by one of the N particles. Particles 
occupying nearest neighbour sites interact with an attractive energy — e < 0. 
The Hamiltonian is then given by 



M p 2 

i=l zm <i,j> 
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where rij = or 1, is the occupation number of the i th site. Mass conservation 
implies 

M 

1=1 

This Hamiltonian can be considered as a simplification of the simple fluid 
Hamiltonian given by Eq. (2). The equivalence of the lattice-gas model and 
the Ising model allows us to use the framework of the latter ( [Lee and Yang| , 



1952| ). An up spin is assigned to an occupied site, a down spin to an empty 



one, and a magnetic field H allows for the conservation (on average) of the 
number of particles. The interaction energy between two nearest neighbour 
up spins is J = e/4. The critical point is located at T = T c , p c = 0.5 in the 
lattice-gas modelj^] and H = in the Ising model. The line defined by H = 
and T < T c in the phase diagram corresponds to the coexistence curve. 

In the Ising model, configurational clusters are sets of nearest neighbour sites 
occupied by an up spin. For a given configuration of spins, the currently called 
"Ising clusters" are determined uniquely, as is shown in Fig. 3. By definition, 
Ising clusters are sets of interacting particles and there is no attraction between 
them. However, it must be emphasized that these clusters have an excluded 
volume which cannot be disregarded at high density. Consequently, unlike the 
Fisher droplets, Ising clusters do not make up a perfect gas of clusters. 




Fig. 3. Ising clusters on a square lattice. There is no (diagonal) bond between the 
next-nearest neighbour spins up. 

The dimension d of configurational space has a crucial effect on thermal and 
geometric phase transitions. Therefore, the verification of the four Fisher cri- 
teria has to be considered according to the dimensionality. 

In two dimensions: Coniglio and his collaborators ( |1977c| ) demonstrate by 
topological arguments that an infinite cluster appears exactly at the critical 
point. Numerical simulations on a square lattice exhibit a power law cluster 
size distribution characterized by an exponent r = 2.1 ±0.1 in agreement with 
Fisher's model (|Stoll et al\ , |1972|) . On the other hand, a renormalization group 

13 By spin up/spin down or particle/hole symmetry, the critical density of the lattice- 
gas model equals to 0.5 regardless of, either the structure of the lattice, or the 
dimension of the configurational space. 
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calculation of exponent v shows that the average radius of the Ising clusters 
diverges like the correlation length, that is v = ( [Klein et ai , 1978 ). The 



first and fourth Fisher criteria are fulfilled by the Ising clusters. 

In contrast, the exponents 7 and (3, which give respectively the critical be- 
haviour of the mean cluster size and the mean size of the largest cluster, do 
not have the expected values of the thermodynamic exponents. Indeed, Sykes 
and Gaunt (|1976| ) show that 7 = 1.91 ± 0.01 (while 7r = 1.75) by means of 
series expansion calculations at low temperature on a triangular lattice. Nu- 
merical simulations confirm this disagreement: 7 = 1.901 ( [Fortuna"to , 2002 ) 



and p = 0.052 ± 0.030 (while (3 T = 0.125) (|Jan et al.\ , [1982|) . In other words 



the second and third Fisher criteria are not verified. In two dimensions Ising 
clusters do not behave like Fisher droplets. 

In three dimensions, none of the four Fisher criteria are met. Along the co- 
existence curve (for H = 0), a percolating cluster appears at a temperature 
lower than the critical one (at a density p < p c = 0.5) ( |(Joniglio| , |1975| ): 



at T ~ 0.94T C and p ~ 0.25 in the cubic lattice from numerical simulations 
( |Muller-Krumbhaar| , |1974a| ,|5|), and at p ~ 0.1 in the face-centered cubic lattice 



from an analytical calculation flSykes and Gaun"t| , |1976|) . 



Let us now consider the behaviour of the Ising clusters in the supercritical 
phase. By using numerical simulations on a square lattice, Odakaki and his 
coworkers ( 1975| ) show that a percolation linep^l divides the supercritical region 



of the T — p phase diagram into two areas: at low density there are only 
small size clusters, compared to the size of the system, and at high density a 
macroscopic percolating cluster connects two opposite sides of the lattice. This 
percolation line is plotted in Fig. 4 for two lattices. In two dimensions, and only 
in two dimensions, the percolation line starts at the critical point. In the limit 
T — > 00, this line tends toward a density equal to the random site percolation 
threshold q c . Indeed, at high temperature, interaction between spins becomes 
negligible compared to thermal agitation. As a result, up spins (particles) 
are randomly distributed with a probability equal to the density and we find 
the particular case of the random site percolation ( Coniglio et all 1977c ). 



As we can see in Fig. 4, the percolation behaviour is almost independent of 
temperature. However, the positive slope shows that a short range attraction 
favors the formation of clusters at low temperature, compared to the random 
site percolation at infinite temperature. 

Critical exponents that characterize the moments of the cluster size distribu- 
tion in the vicinity of the percolation line are equal to the random percolation 
exponents. A calculation using the renormalization group on a triangular lat- 
tice shows that the cluster radius diverges with the exponent v v along the 



1 In a previous work, Stillinger Jr ( 1963| ) had suggested the existence of a percola- 
tion line in the supercritical phase. 
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percolation line ([Klein et al\ , |1978|) . On the other hand, numerical simulations 
provide (3 = 0.145 ± 0.027 (whereas (3 P = 0.139) flOdagaki et al.{ |1975|) . 



T/Tc 




T/Tc 



0.5 



Fig. 4. Schematic phase diagram of the lattice-gas model. The coexistence curve 
(solid line) and the Ising clusters percolation line (dashed line) are represented for 
a square lattice (q c ~ 0.59) on the left, and for a cubic lattice (q c ~ 0.31 ^ on the 
right. For a triangular lattice, q c = 0.5: the percolation line is along the isochore 
p = 0.5 (see text for notation). 

In an infinite Bethe lattice (mean field theory), Coniglio ( |1975| ; |1976| ) demon- 
strates that the percolation line does not start at the critical point, but at a 
lower density than the critical one, like in three- dimension lattices. 



In conclusion, Ising clusters do not behave like Fisher droplets at the ther- 
modynamic critical point, regardless of the dimensionality. Moreover, Ising 
clusters give rise to a percolation line in the supercritical phase. It is impor- 
tant to remark that the presence of a percolating Ising cluster in the super- 
critical phase was considered as non-physical, since on the basis of Fisher's 
model, it would imply thermodynamic singularities ([Reattol, [1970| ; Binder] , 
1975| ; |Muller-Krumbhaar and Stoll| , |1976| ). We shall return to this point later 



when we discuss Swendsen-Wang clusters (see section 4.3). Since Ising clusters 
are too large to behave like Fisher droplets, smaller clusters have to be defined 
(|Binder|, |1976a|). 



4-2 Coniglio- Klein clusters 



Fortuin and Kasteleyn ( |Kasteleyn and Fortuinj , |1969| ; |Fortuin and Kasteleyn 



1972Q propose a correspondence between thermodynamic and geometric quan- 
tities through a formalism based on graph theory within the framework of the 
Ising model. When the probability that two up spins are linked is properly 
chosen, the thermodynamic critical point is exactly a geometric critical point 
characterized by the exponents of the universality class of the Ising model. 

Using an Hamiltonian formulation of random percolation based on an ana- 
lytical extension of the Potts model ( |Wu| , |1982j ) with x colors when x — > 1, 
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Coniglio and Klein ( |1980|) obtain the same results as Fortuin and Kasteleyn, 
and define the so-called "Coniglio-Klein clusters" . For a discussion about the 
relation between the approaches proposed by Fortuin and Kasteleyn and by 
Coniglio and Klein, we refer the reader to ([Hu| , |1984j) and to the Appendix 
of ( |Coniglio| , |2001| ). Two nearest neighbour sites occupied by an up spin are 



connected by a bond with the probability^] 



Pck(Pe) = 1 - e" 



2 



It is sometimes helpful to think of 1 —p c k as the probability of breaking a bond 
between two nearest neighbouring up spins in an Ising cluster. Coniglio-Klein 
clusters are defined as sets of up spins linked by bonds.[^ We remark that the 
probability p c k is independent of the dimensionality. Figure 5 shows a config- 
uration of Coniglio-Klein clusters obtained from the Ising clusters represented 
in Fig. 3. It must be emphasized that unlike Fisher droplets, Coniglio-Klein 
clusters do not form a perfect gas: besides having an hard-core repulsion like 
Ising clusters, they interact through nearest neighbouring up spins that are 
not connected by a bond, as illustrated in Fig. 5. 




Fig. 5. A configuration of Coniglio-Klein clusters on a square lattice obtained from 
the Ising clusters of Fig. 3. The broken bonds between nearest neighbour up spins 
are represented by dashed lines. 

In spite of this crucial difference, Coniglio-Klein clusters satisfy the four Fisher 
criteria at the critical point ( Coniglio and Klein , 1980| ). Numerical simulations 



confirm this result and give for H = and T = T c , values of the geometric 
exponents equal to the thermodynamic ones: in two dimensions, an infinite 
cluster appears at the critical point flUttav^ , |1981| ), 7 = 1.77 ± 0.03 (while 
7T = 1.75), ~ 0.125 = (3 T ( |Jan et a/.| , [1982D , and r = 2.04 (while the Fisher 



exponent is r = 2.07) (|Liverpool and Glotzer , 1996|) . In three dimensions, 
Coniglio-Klein clusters diverge at the critical point (|Stauffer| , |1981| ) with 7 ~ 



15 We choose to express p c k in the lattice-gas framework. Because J = e/4, we have 
in the Ising model: p c k = 1 — e~ 2 ^ J . 

16 Following an idea of Binder ( |1976a| ), Alexandrowicz (1988; 1989) proposes to 
define clusters as sets of up (or down) spins surrounded by a perimeter of up and 
down spins at their average proportion. This leads to the same fractal dimension as 
the Coniglio-Klein 's definition at the critical point. 
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1.2 (while 7t = 1-24), and r ~ 2.25 (while the Fisher exponent is r = 2.21) 
flRoussenq et al\ , |1982| ). 



Coniglio-Klein clusters give rise to a new percolation line, known as the "Kertesz 
line" or "Coniglio-Klein line" , which is shown in Fig. 6. As it should be, this 
line starts at the critical point (p c = 0.5, H = 0) whatever the dimension- 
ality, and reaches the point (p = 1, H = oo) at temperature Tb c given by 

1 — e 2fcTfc c = pb c , where pb c is the random bond percolation threshold. Indeed, 
for p = 1, all the sites are occupied by an up spin (particle) and we find the 
particular case of random bond percolation ( [Heermann and Staufier| , |1981[ ). 




Fig. 6. Phase diagram T — p of the lattice-gas model determined numerically in a 
simple cubic lattice. Coexistence curve (solid line). Critical point (CP.). The Ising 
clusters percolation line (dashed- dotted line) ends at the site percolation threshold 
(q c ~ 0.31) when T —* oo. The Coniglio-Klein percolation line, or Kertesz line 
(dotted line), ends at the bond percolation threshold when p = 1 (from Campi arid] 
Krivinj \l9^ ). 

As we can see in Fig. 6, the supercritical area of the phase diagram is divided 
into a non-percolating region at high temperature and low density, and a 
percolating region at low temperature and high density. Along the Kertesz 
line, geometric exponents are equal to those of random percolation ( |Coniglio| 



and Klem| , [L980| ; (Staufier| , |1981| ). Therefore, values of geometric exponents 
pass disco ntinuously from the thermodynamic ones at the critical point, to 
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the random percolation values along the line, for T > T c and p > p c (H > 0). 



It is interesting to remark that when interactions with next or second neigh- 
bours are taken into account, the critical point remains a geometric critical 
point for Coniglio-Klein clusters^ 7 ] ( |Jan et aT\ , |1982j ) . This result suggests that 
moving on to the continuum limit and considering a more realistic potential 
with a larger range of attraction, the critical point is a percolation thresh- 
old.^] This possibility, as well as the existence of a percolation line in a real 
fluid, will be discuss later on in section 5. 



To close this subsection, we briefly discuss the behaviour of Coniglio-Klein and 
Ising clusters below the coexistence curve. At thermal equilibrium, Monte- 
Carlo simulations in three dimensions show that no critical percolation be- 
haviour exists when phase separation occurs: critical percolation lines do not 
get into the two-phase region ( |Campi et al\ |2002a| ). Nevertheless, below the 
coexistence curve, there is always a macroscopic cluster (either a Coniglio- 
Klein or an Ising cluster), in the sense that a finite fraction of the particles 
belongs to the largest cluster.^] This is not the purpose of this review ar- 
ticle to deal with the out of equilibrium and nucleation properties of simple 
fluids, however, we mention that Coniglio-Klein clusters, as well as Ising clus- 
ters, fulfill the phenomenological laws of Classical Nucleation Theory below 
the coexistence curve flHccrmann et at , 1984 ; [Ray and Wang , 1990 ). On the 
other hand, in mean field theories, the compressibility diverges not only at 
the critical point but also along the well-defined spinodal line which separates 
the metastable and unstable regions. The mapping between thermodynamics 
and geometry in terms of Coniglio-Klein clusters cannot be extended into the 
metastable region flSchioppa et al.j |1998|) . However, by introducing the bond 
probability Pg = 1 — e~^ l ~ p ^ between particles, the mean size of these new 
clusters diverges along the mean field spinodal line ( Heermann et af\ 1984 ; 
[Klein], HOT. 



17 For Ising clusters, the density of the geometric critical point along the coexistence 

|1985| ). 



5ug et al. 



curve obviously decreases as the range of attraction increases 

18 As a matter of fact, Fortuin and Kasteleyn's formalism has been generalized to 
the case of a short range potential in a continuous space ( prory , 1996| ). 

19 To put it more precisely, the macroscopic cluster can be either too compact to 
span the system or enough spread out to percolate, according to temperature and 
density values. In the last case, although the largest cluster connects two opposite 
sizes of the lattice, we recall that no critical percolation transition is observed, the 
distribution of finite size clusters remains exponentially decreasing and does not 
display any power law behaviour flCampi et al\ , 2002a ). 
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4-3 Swendsen-Wang clusters 



As we have seen in the last section, Coniglio-Klein clusters behave like Fisher 
droplets at the critical point (T = T c , H = 0). However, the percolation line 
in the supercritical phase (T > T c , H > 0) was interpreted as a non-physical 
manifestation of these clusters. Within the context of the Fisher droplet model 
(or more generally in the perfect gas of clusters model), a percolation be- 
haviour should imply singularities of thermodynamic quantities, which is for- 
bidden in the supercritical phase by the theorems of Yang and Lee ( |1952|) . The 
aim was then to define clusters that have a geometric critical behaviour (per- 
colation threshold) at, and only at, the thermodynamic critical point. In this 
spirit, the Kertesz line was rather cumbersome ( |Coniglio et al.\ , |1989| ; Kerteszj , 
1989| ; jjtautteil , |I990| ). 



/ / 
/ / 
/ / 
/ / 





Fig. 7. Up Spins (particles) are represented in the lattice- gas formulation by a black 
dot. The white dot is the "ghost spin". Bonds between up spins are active with 
probability p c k = 1 — e~ 2 @ J (heavy line) and bonds between up spins and the "ghost 



spin" are active with probability ph = 1 - 
connected only through the "ghost spin". 



-2/3H 



(dashed line). Some up spins are 



In order to get rid of the percolation line, Coniglio and his collaborators ( 1989|) , 
as well as Wang (|1989|) , introduce a definition of clusters based on the work of 
Fortuin and Kasteleyn QKasteleyn and Fortuinj , |1969| ; |Fortuin and Kasteleyn , 
1972j ), and Swendsen and Wangf 20 ] ( |1987| ). These so-called "Swendsen-Wang 
clusters" are defined like the Coniglio-Klein clusters, but in the presence of 
a positive (negative) magnetic field, H, up (down) spins are in addition con- 
nected to a "ghost spin" with the probability ph = 1 — e - '' 1 ', where h = 2(3H. 
This model is best explained by referring to Fig. 7. At each point of the phase 
diagram, the cluster finite size distributions of Swendsen-Wang clusters nf", 



20 Swendsen-Wang algorithm, which allows one to reduce the critical slowing down 
at the critical point, is essentially based on the Fortuin and Kasteleyn's formalism 
(|Swendsen and Wang , 1987 ; Wang and Swendsen , 1990 ; [Liverpool and Glotzei , 
119961). 
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and Coniglio-Klein cluster n c s k , are related by ( |Wang| , pL989|) : 



n 



sw -= 2n c S~ Ws where s is finite. 



Hence, as soon as H is different from zero, a percolating cluster, which con- 
tains the ghost site, is present in the system whatever the temperature. In- 
deed, for H > 0, even if the magnetic field is weak, a macroscopic fraction 
(1 — e~ 2/3H ~ 2/3H) of up spins is connected to the "ghost spin". As a re- 
sult, there is a percolating cluster at each point of the phase diagram located 
above the coexistence curve ( [Stauffer' , 1990 ; Adler and Stauffcr , 1991| ). In the 



thermodynamic limit, any cluster that contains the "ghost spin" is infinite, 
and vice versa. The infinite cluster is unique. Finite Swendsen-Wang clusters 
behave as we expected: they verify the four Fisher criteria at the critical point, 
as the Coniglio-Klein clusters do. Besides, their mean size does not diverge in 
the supercritical phase. The percolation line is eliminated. 

However, even if the finite Swendsen-Wang cluster size distribution does not 
have a power law behaviour in the supercritical phase, something happens 
along the Kertesz line. Wang ( |1989|) calculated nf", for large but finite clus- 
ters, by means of numerical simulations. According to the region of the phase 
diagram (see Fig. 6), he found: 

• Above the Kertesz line, n s s w ~ e -cis-\h\s 

• Along the Kertesz line, n s s w ~ s~ T e~' ft ' s 

2 

• Below the Kertesz line, nV" ~ e - c 2 s ^-\ h \ s 



where the coefficients c\ and C2 depend on temperature and magnetic field. The 
Swendsen-Wang cluster size distribution decreases exponentially everywhere 

2 

in the supercritical phase. However, the "surface term" of n s s w , varying like ss 
below the line, goes to zero along and above the line. Kertesz ( |1989j ) and Stauf- 
fer ( |1990| ) suggested that the cancellation of the surface term in the expression 
for the cluster size distribution should imply a weak singularity of the cluster 
free energy.^] Discussing this result, Stauffer and Aharony ( [1994]) wrote: "If 



correct, it would mean that we have to make more precise the century-old 
wisdom that nothing happens if we move continuously from the vapor to the 
liquid phase of a fluid by heating it above the critical temperature." 

At the beginning of the nineties, this claim sounded conclusive. Nevertheless, 
we believe this interpretation of the Kertesz line to be incorrect and propose a 
new point of view. Indeed, it must be emphasized that this argument is based 



21 According to Kertesz ( 1989 ) the contribution of the percolating cluster to the 
total free energy should cancel out this weak singularity and the partition function 
would be of course analytical. 
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on the hypotheses of the perfect gas of clusters in which (see Eq. (8)) 

n s = q s z s = e ~^ +slnz (19) 



Remember that Coniglio-Klein clusters, as well as Swendsen-Wang clusters, 
do not form a perfect gas of clusters. We can see therefore no valid reason 
for writing their cluster size distributions as Eq. (19). In other words, the 
particular behaviour of the Coniglio-Klein and Swendsen-Wang cluster size 
distributions along the Kertesz line does not imply even weak thermodynamic 
singularities. Furthermore, as far as we know there is no physical interpretation 
of the ghost spin and Swendsen-Wang clusters. We shall see that this is not 
the case for Coniglio-Klein clusters. 

Before we move on to the next section, we point out a general fact about 
clustering in simple fluids: the definition of clusters does depend on the phys- 
ical problem we want to study. For example, Ising clusters do not provide a 
geometric interpretation of the liquid-gas transition. However, properties of 
electrical conductivity can be understood in terms of these clusters of inter- 
acting particles ( |Odagaki et al.\ , |1975| ; |Yang and Thompson| , |1998| ). Moreover, 



in the continuum case, configurational clusters have been used to describe the 
dynamics of hydrogen bonds in liquid water ( [Luzar and Chandler| , |1996| ; |Stari 
~Wal\ [T9991) . 



5 Self-bound clusters 



From now on, our purpose is to define clusters as self-bound sets of particles. 
A typical physical situation in which it is relevant to introduce "self-bound 
clusters" is the fragmentation of a piece of matter (such as an atomic nucleus, 
an atomic aggregate, or a liquid droplet) ( Campi et al.\ , [2001aj ). In particular, 
a question of interest is the existence of physical clusters in a dense medium, 
before fragmentation occurs. This topic seems to contrast strongly with the 
discussion of the last sections. Yet, we shall see how such clusters can shed 
some light on the geometric interpretation of the liquid-gas phase transition. 

Self-bound clusters would be stable if they were isolated from each other. In 
a dense medium, their definition cannot rest only on a criterion of proxim- 
ity in configurational space. Such a definition would make no allowance for 
the relative motion of the particles. At high temperature, two particles that 
strongly interact are close to each other at a given time. Because of their high 
velocities, they may quickly move away from each other. A definition of clus- 
ters that takes into account the relative velocities of the particles corrects this 
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inadequacy.^] 



In section 5.1, we present two definitions based on energetic criteria which 
allow us to identify self-bound clusters in a fluid knowing the positions and 
velocities of the particles. The first definition is global and based on minimiz- 
ing the interaction energy between clusters. The second one is a local energetic 
criterion of bond activation. Next, we show that Coniglio-Klein clusters are 
self-bound on average. This crucial link between thermodynamics and the 
morphology of a fluid in the lattice-gas model will provide a geometric inter- 
pretation of the liquid-gas phase transition in terms of physical clusters. This 
point will be developed in section 5.2 for a Lennard- Jones fluid. 



5.1 Definitions of self-bound clusters 



5.1.1 A global criterion 

At very low density, clusters are naturally defined as sets of particles that 
strongly interact with each other and do not interact (or only very loosely) 
with particles belonging to other clusters. This is the domain of validity of the 
perfect gas of clusters model. 

In a dense fluid, by an argument of continuity, self-bound clusters can be 
defined by the partition of particles that minimizes the interaction energy (in 
absolute value) between clusters, or equivalently that minimizes the sum of 
the clusters internal energies. Indeed, the total energy E of the system divided 
in an arbitrary partition of clusters can be written 

E = J2 \mN c V% + Y^Uc+J2 V UC, C) (20) 
c 1 c c,c 



where J2c * s a sum over the clusters, and Nc, Uc, and Vc are respectively 
the number of particles, the internal energy, and the velocity of the center 
of mass of cluster C. The term Vi n t(C, C) is the interaction energy between 
clusters C and C . We then have to find the partition of clusters that minimizes 

U(P)=Y,U C (21) 

c 



This minimization problem has an exact, but not necessarily unique, solution. 



22 An alternative approach is to introduce an arbitrary bond life time, in addition to 
the cutoff distance: two particles are linked if they remain within a certain distance 
during a given time (Pugnaloni et at, 2000| ). However, not one but two adjustable 
parameters have to be chosen. 
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When the number of particles is small enough to enumerate all the possible 
partitions, a solution can be found easily. Otherwise, we can use two different 
methods for finding an optimal solution. 

The first one was proposed by Dorso and Randrup (|1993| ). It is based on 
an optimization procedure called "Simulated Annealing" (|Kirkpatrick| , |1984|) : 
consider an arbitrary partition of clusters P, a particle is randomly chosen 
and transferred from its original cluster C to another cluster C . Of course, 
the particle remains at the same position in the fluid but we decide that it 
now belongs to another cluster. We obtain a new partition P . The variation 
of the total internal energy between the two partitions is 

AU = U(P') - U(P) = U' C + U' c , -U c - U c > 

where U' c and U' c , are respectively the internal energies of clusters C and C 
after the particle has been transferred. If AU < 0, the partition P is accepted. 

' At/ 

If AU > 0, the partition P is accepted with the probability e T ^ , where T e is 
a control parameter. The value of T e is progressively decreased until there is no 
more modification of the cluster partition on a given number of attempts. As 
usual with Metropolis methods, the total internal energy of the final partition 
may be only a local minimum. 

A second method, denoted "BFM" (for Binary Fusion Method) has been de- 
veloped by Puente (|1999|) : initially, we suppose that there are only monomers 
in the fluid. The multiplicity, which is the number of clusters, is maximal, that 
is mo = N. Then, at each iteration we merge the two clusters that provide the 
lowest value of the function U(P) (see Eq. (21)). The multiplicity is then re- 
duced by one at each iteration. Therefore, at the first iteration a dimer is built, 
at the second either another dimer or a trimer, and so on. The process stops 
when the total internal energy cannot be reduced any more. This determinis- 
tic method follows a certain trajectory in the phase space of the partitions of 
clusters that does not necessarily lead to the optimal partition. On the other 
hand, the Simulated Annealing method can probe the phase space in a prob- 
abilistic way, which depends on arbitrary chosen parameters (like the number 
of Metropolis trials or the way T e decreases). According to their authors, these 
two procedures provide clusters that are stable against particle emission. The 
definition of stability against particle emission is presented at the beginning 
of appendix B. 

By means of numerical simulations of a Lennard- Jones fluid, Puente ( p.999|) 
has shown that the Simulated Annealing and the BFM methods produce very 
similar cluster size distributions, but the last one is 10-20 times faster for the 
system considered. However, it has to be noted that both of these methods are 
very time consuming and cannot be used in practice for large systems made 
up of more than (say) N = 1000 particles. 
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5.1.2 A local criterion: Hill clusters 



The definition of clusters proposed by Hill (|1955|) rests on a bond activation 



between pairs of particles. At a given time, two particles i and j, of respective 
velocities Vi and Vj, are linked if their relative kinetic energy K r is lower than 
the negative of their attractive interaction energy u(rij): 



K r = -^{vi — Vj) 2 < —u(rij) i and j are linked. (22) 



In other words, there is a bond between two particles if they are close enough 
to each other in the phase space, and not only in the configurational space. 
A Hill cluster is by definition a set of particles connected by bonds two by 
two. When the positions and velocities of the particles are known, for instance 
by means of molecular dynamics simulations ( [Allen and Tildesleyj , |1987| ), this 
local definition is completely deterministic and does not depend on adjustable 
parameters. It has to be noticed that the Hill's criterion can be used as well 
with Monte-Carlo calculations, although the velocities of the particles are not 
provided. Indeed, we checked by means of molecular dynamics simulations 
that mean cluster size distributions are not modified when we replace the 
velocity of each particle, provided by molecular dynamics simulations, with 
velocity components taken at random from a Gaussian distribution character- 
ized by the temperature of the system (|Satorj |2000|) . As far as mean cluster 



size distributions are concerned, position-velocity correlations are irrelevant 
and Monte-Carlo calculations can be used by allocating to each particle a 
velocity from a Gaussian distribution. Of course, time dependent properties 
and dynamics of the self-bound clusters are investigated only by molecular 
dynamics simulations. 

A question that naturally arises at this point is the comparison between the 
global (Simulated Annealing and BFM) and the local (Hill) criterion. Perform- 
ing molecular dynamics calculations with a Lennard- Jones fluid, Campi and 
his collaborators ( Campi et ai , 2001a ; Sat or , 2000 ) have found that the clus- 



ter size distributions obtained with these three definitions are very similar. To 
quote Hill (|1955|) : "A cluster is [. . .] a fairly definite physical concept so that 



all physically reasonable definitions should lead to rather similar predictions 
of the cluster size distribution." However, for a system of N = 1000 particles, 
Hill's criterion is 300 times faster to implement and can be used in practice 
for much larger systems, as we shall see in section 5.2. 

Because it is arduous to carry out an analytical approach by using the Hill's 
criterion, it may be relevant to deal with a probabilistic formulation of the 
criterion (22). That requires us to calculate the probability that two particles 
are linked. 
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In a fluid at equilibrium, velocity components of the particles are given by 
a Gaussian distribution characterized by a variance l/(/3m), where m is the 
mass of the particles and (3 = 1/kT. The criterion (22) involves the relative 
kinetic energy between particles i and j which can be seen as the kinetic 
energy of a fictitious particle of reduced mass fi = m/2 and relative velocity 
v r = vi — Vj. The velocity components of the fictitious particle are then given 
by a Gaussian distribution with variance = l/f3m. 

The probability pmuij", (3) that two particles separated by a distance r are 
linked is then the probability that K r = < —u(r), that is 



r y/-Pu(r) 2 2 , 
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where T(n,x) is the incomplete gamma function. In this probabilistic frame- 
work, bonds are independent of each other, like in a random percolation pro- 
cess. A a consequence, we lose the velocity correlations between particles which 
are central in the deterministic criterion. We shall see in section 5.2 that clus- 
ters defined with the deterministic definition are much larger than the proba- 
bilistic clusters. 

By using the probability PHiii( r , P), Hill ( |1955[ ) introduced two effective poten- 
tials between bound and unbound pairs of particles for expressing the mean 
number of clusters as density and activity series expansions.!" 2 ^ A very similar 
approach has been used recently to calculate the first order corrections (namely 
by considering the interaction between a monomer and a self-bound cluster) 
to the perfect gas of clusters at low density ( [Soto and Corderc] , |1997| ) . Coniglio 



and his collaborators (|1977a| ; |1977b| ) extended the Hill's method and devel- 
oped an integral equation theory of the pair connectedness function, which 
is proportional to the probability that two particles belong to the same clus- 
ter f^j With an Ornstein-Zernike-type equation and a closure relation, like the 
Percus-Yevick approximation, the mean cluster size can be obtained from the 
integration of the pair connectedness function. Very recently, a generalization 
to the case of the deterministic Hill's criterion has been proposed and should 
provide interesting results ( Pugnaloni et all |2000| . |2002|) . 



23 This method, which takes into account the interactions between clusters, can be 
formally applied to any microscopic definition of clusters. 
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For a comprehensive review about integral equation theory, see (|Caccamo , 1996 ). 
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5.1.3 Coniglio- Klein and self-bound clusters 

In the lattice-gas model, interaction is equal to — e between nearest neighbour 
particles, and otherwise. The probability that two nearest neighbour particles 
are linked becomes flCampi and Krivinej , |1997| ) 



4 r 2 
Pmu(Pe) = —= / e~ x x 2 dx 
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(23) 



We show in appendix B that probabilistic Hill clusters are stable on average 
against particle emission. By performing a Laguerre expansion of Eq. (23), 
Campi and Krivine ( 1997| ) have found 



Pmu(pe) = 1 - 0.911e-r - 0.177(1 - (3e)e~^ + ... 
= Pck(M - 0.088(1 - 2.01/3e)e - ^ + . . . 

^Pck((3e) (24) 

where p c k(f3e) = 1 — e~^~ is the probability introduced by Coniglio and Klein 
(see section 4.2). Equation (24) shows that probabilities p c k and pmu are very 
close to each other (for example at the critical point of the cubic lattice, 
the relative error between the two probabilities is less than 6%) J 25 ] This is a 
crucial result: Coniglio-Klein clusters, which have been introduced specifically 
to provide a correspondence between geometric and thermodynamic critical 
behaviour, do have a physical interpretation in terms of self-bound clusters. 
Conversely, self-bound clusters give rise to a percolation line which starts at 
the critical point in the framework of the lattice-gas model. 



5.2 Percolation lines in simple fluids 



Numerical simulations of a Lennard- Jones fluid provide a more realistic theo- 
retical framework than the lattice-gas model. With the help of canonical Monte 
Carlo and microcanonical molecular dynamics calculations ( |Allen and Tildes^ 
[ey| , |1987|) , Campi and his collaborators dSatoij [2000| ; |Campi et a/.| , pOOlaf) 



have studied systems made up of a large number of particles [N = 11664) 



25 Pan and Das Gupta ( |1995 ) had already noticed this similarity using numerical 
simulations. 
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interacting through a Lennard- Jones potential given by Eq. (1). The phase 
diagram of this simple fluid is presented in Fig. 8. It must be emphasized that 
the phase diagram, as well as the results we present below, are qualitatively in- 
sensitive to this particular choice of potential and can be considered as generic 
for simple fluids. To reduce the computer time calculation, self-bound clusters 
are recognized by using the local deterministic criterion proposed by Hill (see 
section 5.1.2). 
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Fig. 8. Phase diagram of the Lennard-, J ones fluid. The dashed line indicates the 
coexistence curves between the gas (G), liquid (L) and solid phases (S). The contin- 
uous line is the critical percolation line of self-bound clusters. Temperature T and 
density p are in units of e and a~ 3 respectively (see Eq. (1)) (from Campi et al 



2001a) 



As a first result, a "macroscopic" cluster appears as soon as the liquid-gas 
phase transition takes place. Figure 9 represents two cluster size distributions 
plotted around the point (e) of the phase diagram (see Fig. 8), just above 
and just below the coexistence curve. We can see a drastic change in the 
mean number of large clusters which corresponds very well to the crossing 
of the coexistence curve. Obtained without any free parameter, this clean 
correspondence between thermodynamic and morphological changes cannot 
but support the physical nature of these clusters. It is important to note that 
along the coexistence curve, the cluster size distribution is not a power law 
and therefore, does not present a geometric (percolation) critical behaviour. 

What is more, self-bound clusters give birth to a critical percolation linep^l 



26 This line should not be confused with the Fisher- Widom line ( Fisher and Widonj 
|1969 ), which divides the supercritical region into a gas-like and a liquid-like domain 
(with probably no relevance to clustering) according to the behaviour of the radial 
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Fig. 9. Cluster size distributions n s at points (a)-(b), (c)-(d) and (e) of the phase 
diagram (see Fig. 8). For curve (d): the contribution of the percolating cluster is 
sharply peaked around s ~ 9500. For curves (e): T con d is the temperature of conden- 
sation at p = 0.1 (from Campi et ai, 200 lap . 



represented in Fig. 8. To illustrate this percolation behaviour, the cluster size 
distribution is plotted in Fig. 9 for several points of the phase diagram. Along 
the line the cluster size distribution behaves as a power law characterized by 
an exponent r = 2.20 ± 0.05 (while in random percolation theory r p = 2.18). 
The fractal dimension of the clusters, as well as the other exponents that 
characterize the moments of the cluster size distribution in the vicinity of 



distribution function, nor with the extrapolation of the rectilinear diameter line 
( Nishikawa and Morita| , 2000] ), which is the locus of points of maximum compress- 
ibility. 
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the line, are also in good agreement with those of random percolation theory 
(ICampi et alj gjfflg ; [Satorl g000|) . Yet, this result was found on a deterministic 



dynamical framework, without any explicit reference to a random (site or 
bond) percolation mechanism. This is even more striking when the global 
definition is used instead of the bond activation criterion. 

The percolation line starts at the critical point, within error bars due to the 
finite size of the system and to critical slowing down. Just at this point, for 
the same technical reasons, geometric critical exponents were not evaluated 
with enough accuracy to distinguish between the universality classes of the 
Ising model and of the random percolation. Furthermore, the total energy of 
the system remains almost constant along the percolation linef 27 ] ( pampi ei 



a I. 



], |2001a| ). The origin of this energy invariance, which results from a subtle 



balance between internal, center of mass kinetic, and inter-cluster potential 
energies (see Eq. (20)), is not understood flCampi et al.[ [2001b|) . 



It is interesting to notice that along the condensation curve, and in particular 
at the critical point, self-bound clusters seem to behave like Fisher droplets. 
Hence, their cluster size distributions can be fitted by Eq. (13), which was 
proposed by Fisher. However, self-bound clusters do not form a perfect gas 
as Fisher droplets do, and there is absolutely no reason for writing their free 
energy as —f3F s = In n s — s In z. Indeed, the contribution of the interaction be- 
tween self-bound clusters to the total potential energy of the system is at least 
40% along both the percolation line and the condensation curve ( |Sator| , |2000| ). 



Even at very low density, such as p ~ 0.05 in two dimensions, the interaction 
between a single monomer and a self-bound cluster markedly modifies the clus- 
ter size distribution of the ideal gas of clusters ( [Soto and Corderol , |1997|) . In 



addition, we recall that the definition of the self-bound clusters tends to mini- 
mize interaction between clusters. These results show strikingly that whatever 
microscopic definition of clusters is chosen, it is not realistic to treat them as 
a perfect gas of clusters, even at low densities ~ p = p c /3 ~ 0.1. 

For a given cluster size s, the degrees of freedom associated with the motion of 
the particles inside the clusters, and the center of mass motion of the clusters, 
are respectively characterized^] by an "internal effective temperature" T*(s) 
and a "translational effective temperature" T cm (s). Analytical calculations 
(|Soto and (Jordero| , |1998a| ,[E|) as well as numerical simulations flCampi et ai\ , 



2001b| , |2002b|) show that whatever the size of the clusters and the thermody- 



namical state of the system, T*(s) < T < T cm (s) where T is the real thermo- 
dynamic temperature of the system. Simply put, the system can be seen as a 



27 The same observation has been made about the Kertesz line in the lattice-gas 
model ( pampi et al , 1999| ). 



28 These effectives temperatures, T*(s) and T cm (s), are defined as twice the corre- 
sponding average kinetic energy divided by the number of degrees of freedom, that 
is 3(s — 1) and 3 respectively in three dimensions. 
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"hot" gas of "cold" physical clusters. This result, due to the self-bound nature 
of the clusters, is particularly relevant when studying the fragmentation and 
expansion of a small system. Indeed, once the clusters are isolated from each 
others, the internal effective temperature becomes the real temperature of a 
fragment in the microcanonical ensemble ( |Campi et al.\ , [2002 b| ). 



As we have seen in section 5.1.2, self-bound clusters can also be defined by a 
simplified probabilistic criterion. By means of molecular dynamics simulations 
of a Lennard- Jones fluid, Pugnaloni and his collaborators ( |2002|) have recently 
shown that this definition strongly overestimates the percolation density by 
a factor 1.5-2, especially at high temperatures. In other words, probabilistic 
clusters are much smaller on average than the deterministic Hill clusters. On 
the other hand, the mean cluster size of the probabilistic Hill clusters can be 
calculated using an integral equation theory^] (|Coniglio et al\ |1977a|) . The 
percolation density is then determined by the divergence of the mean clus- 
ter size. When compared to the molecular dynamics results, the percolation 
line of the probabilistic clusters is shifted to much lower densities and may 
seem to roughly coincide with the percolation line of the deterministic clusters 
flConiglio et g] , |H)77a| ; |Pugnaloni efal\ PUTI^ ). 




Fig. 10. Percolation lines of the configuration^ clusters in the T—p phase diagram of 
a Lennard- Jones fluid determined by molecular dynamics simulations for two system 
sizes: N = 108 (black points) and N = 500 (circled points). The corresponding value 
of the cutoff distance b is indicated for each percolation line. Note that b m i n — 1.12 
is the distance corresponding to the minimum of the Lennard- Jones potential (from 
Heyes and Melro~se\ \l98i ). 



For completeness, we mention that configurational clusters (see section 4), 
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Xu and Stell ( 1988| ) used this approach to study the percolation behaviour of 



a Yukawa fluid. However, clusters are defined either by a cutoff distance or by a 
probability of connection which depends on an adjustable parameter. 
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defined by a cutoff distance b, also generate a percolation line which almost 
follows an isochore (p = cte), as do Ising clusters in the lattice-gas model 
flGeiger and Stanley] , |1982| ; |Bug et aJ\ |1985| ; |Safran et al\ , |1985| ; [Xu and Stell 



|1988| ). Of course, the position of the line in the phase diagram depends cru- 
cially on the choice of b flHeyes and Melrose) , gggg , gggg |Xu and SteTj , pJ88|) , 
as can be seen in figure 10. 

In the very particular case of Baxter's sticky-hard-sphere model (or adhesive- 
hard-sphere model) (|Baxter| , |1968|) , configurational clusters become identical 
to self-bound clusters. The potential is a hard spheres plus a square-well po- 
tential, in the limit where the depth of the well goes to infinity, and its width 
goes to zero, in such a way that the product of depth and width remains con- 
stant. The cutoff distance b is then obviously the diameter of the core. Besides, 
once in contact, two particles are permanently bound. Baxter ( pL968| ) solved 
this model analytically, using the Ornstein-Zernike equation in the Percus- 
Yevick approximation, and showed that the system undergoes a liquid-gas 
phase transition (see Fig. 11). On the other hand, Chiew and Glandt ( j!983|) 
used the approach proposed by Coniglio and his collaborators ( |1977a| ) (see 
section 5.1.2) to locate the percolation line in the phase diagram of the sticky- 
hard-sphere model. As is shown in Fig. 11, this line does not pass by the critical 
point. However, following Kranendonk and Frenkel (|1988|), we note that its 



location for extreme values of the volume fraction^] casts some doubt upon 
the adequacy of the Percus-Yevick approximation. Indeed, the percolation line 
starts at = 0, predicting the existence of a percolating cluster in vacuum, 
and goes to the unphysical value 0=1, greater than the close packing vol- 
ume fraction. In contrast, Monte-Carlo simulation predictions show that the 
percolation line starts very close to the critical point and goes through the su- 
percritical phase as in the case of a Lennard- Jones fluid flSeaton and Glandt| , 
|1987a| ; [Kranendonk and Frenkel , |1988j ) . According to |Kranendonk and Frenkel 
( |1988| ), the prolongation of the percolation line under the coexistence curve 



should be taken with caution. 



6 Summary and outlook 



Ever since Mayer proposed his theory of condensation in 1937, the correspon- 
dence between thermodynamics and the morphology of simple fluids in terms 
of clusters has attracted a lot of interest. In this last section, we summarize 
the various stages of this research area, and close by mentioning some open 
questions. 



30 The volume fraction is related to the density by = irpa 3 /6, where a is the 
diameter of the particles. 
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Fig. 11. Phase diagram of the sticky-hard- sphere model. The stickiness parameter r 
(vertical axis) is a measure of the temperature and <p is the volume fraction (hori- 
zontal axis). The dashed line represents the spinodal line in the Percus-Yevick ap- 
proximation. The percolation line determined by the Percus-Yevick approximation, 
dotted line ( Chiew and Gland\ , \1983j) , and by numerical simulations, black squares 
(Kranendonk and Frenkei, 1988), and crosses ^Seaton and Gland\ , 1987a\ j^). The 
horizontal bars give the positions of percolation thresholds below the spinodal line 
(from Kranendonk and Frenkei, 1988\). 



Initially, the purpose was to introduce clusters as tools in order to calculate the 
thermodynamic quantities of a classical fluid, and in particular to study the 
liquid-gas phase transition and the critical point. To this end, condensation is 
defined as the formation of a macroscopic cluster. By considering a real fluid 
as a gas of non-interacting clusters, one can write the equation of state of the 
system as a function of the cluster size distribution. Consequently, the mean 
cluster size varies like the compressibility and then should diverge at, and 
only at, the thermodynamic critical point. This is the case by construction for 
Fisher droplets which provide a morphological description of the critical point. 
But these phenomenological non-interacting clusters cannot be identified as 
sets of particles in the fluid. The next step was to find a microscopic definition 
of clusters that behave like Fisher droplets. 
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In the lattice-gas model, the problem amounts to finding a definition of clus- 
ters that produces a percolation threshold at the critical point. Ising clusters, 
defined by the simple criterion of proximity in real space, do not fulfill this 
condition. In contrast, Coniglio-Klein clusters, whose definition is based on 
the mathematical work of Fortuin and Kasteleyn, are explicitly designed to 
this end. But this is not the end of the story: Coniglio-Klein clusters exhibit 
a percolation behaviour not only at the critical point, but also along the so- 
called "Kertesz line" which goes through the supercritical region of the phase 
diagram. Within the context of the perfect gas of clusters model, this per- 
colation line should imply forbidden thermodynamic singularities. However, 
because Coniglio-Klein clusters do not form a perfect gas of clusters, there 
is no valid reason for eliminating the Kertesz line. It remains to give a phys- 
ical interpretation of this percolation line. This is what Campi and Krivine 
have done by showing that Coniglio-Klein clusters are self-bound on average. 
This interpretation in terms of self-bound clusters opens up a new field of 
investigation. 

Self-bound clusters can be studied in the more realistic framework of a Lennard- 
Jones fluid. It is interesting to notice that these clusters were introduced not to 
provide a correspondence between thermodynamics and geometry, but rather 
to describe processes like the fragmentation of a piece of matter. Neverthe- 
less, self-bound clusters strikingly present the morphological manifestations 
associated with the liquid-gas phase transition, and therefore can be consid- 
ered as the microscopic physical counterparts of the phenomenological Fisher 
droplets. What is more, self-bound clusters give rise to a percolation line in 
the supercritical phase. An important open question, which arises from their 
definition based on physical grounds, is whether this line could be observed in 
real fluids. 



Various real systems behaving like simple fluids can be used to search for 
signals of the self-bound clusters and the percolation line. For example, the 
effusion of a fluid through a pin hole allows one to infer the cluster size distri- 
bution from the mass yield of escaped clusters with a given velocity. Hence, 
the presence of clusters (dimers and trimers) was found in the gas phase at 
very low density ([Miller and Kusch] , |1955| , |1956|) . However, this kind of ex- 
periment seems to be inappropriate to explore a fluid at high density, in the 
region of the percolation line. On the other hand, experiments dealing with 
binary fluids (which belong to the universality class of the 3d Ising model 
( |Chen et at.\ |1983| )) were performed to give a morphological description of 
the critical point in terms of clusters ( puenoun et al\ , |1989|) . The relationship 
between these clusters, defined from persistent density fluctuations and the 
self-bound clusters is not yet understood. Very recently, a sharp decrease of 
H + ion mobility in a mixture of isobutyric acid and water has been ascribed 
to a percolation line of "dynamic clusters" which starts at the critical point 
( |Bonn et al\ |2002[ ). At first glance, the vertical position of this curve in the 
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phase diagram would suggest a closer connection with configurational clusters 
(see Fig. 5 in ([Bonn et al , 20021 )) . However, the slope of the percolation line 
of the self-bound clusters may depend on the system considered. The question 
of whether these "dynamic clusters" are self-bound, still awaits an answer. 

From a rather different perspective, the relation between clustering and ther- 
modynamics is critical to the determination of the phase diagram of small 
systems, like atomic nuclei and aggregates ( pill , |1963| ). In a first approxi- 
mation, the nuclear interaction can be seen as a two body potential with a 
repulsive hard-core and a short-ranged attraction. By disregarding quantum 
effects at high enough energy, hot nuclear matter is thought to behave like a 
classical simple fluid QBondorf et al. 



1995|; |Das Gupta et g[\ |2001|, and ref- 



erences therein). A crucial question which then arises is whether traces of a 
liquid-gas phase transition can be observed in an atomic nucleus ([Siemens] , 
|1983j ) ? A direct observation of such a signal is of course quite arduous in a 
small system, and today there is still no conclusive answer to this question. 

To explore the phase diagram of nuclear matter, atomic nuclei are heated and 
compressed by collisions at high bombarding energy (about 100 MeV/nucleon) 
with other nuclei. As a result, an excited nucleus expands and breaks into 
fragments which can be detected using 4tt multidetectors (see references in 
Bondorf et al.[ |1995|) . Basically, only the fragment size and the kinetic energy 
distributions can be directly measured. The challenge is then to infer the 
thermodynamical state of the system before fragmentation occurs from these 
two distributions which characterize the asymptotic fragments. 



To this end, several theoretical approaches have been proposed. "Statistical 
Equilibrium Models", like the SMM ([Bondorf et~al\ p95[) and the MMMC 
(|Gross| , |1997|) models, assume that at a given stage of the expansion the sys- 
tem is an ensemble of non-interacting spherical fragments. Using the perfect 
gas of clusters model, and some elements of nuclear physics, these models 
have been successful in describing the observed fragment size distributions 
QBotvina et al.\ , |1995| ). However, the kinetic energies are too low compared to 
experimental data, and worst of all, the strongly restrictive hypotheses of the 
perfect gas of clusters prevent exploration of the high density region of the 
phase diagram. On the other hand, classical molecular dynamics simulations 
of the expansion and fragmentation of a small system allow one to investigate 
the whole phase diagram ( [Schlagel and Pandharipandej , |1987| ; [Belkacem et al. 



1995| ; |Chernomoretz et al.\ , |2001| ; |Campi et al.\ , [2002 b| ). In this way, it has been 



shown that self-bound clusters in the hot and dense fluid are the precursors 
of the observed fragments ( Dor so and Kandrup , 199cl ; Campi et al. , 200 lb| . 
2002b|) . In particular, self-bound clusters in the initial system and asymptotic 
fragments have the same size distribution. The initial thermodynamic state 
could be deduced from the measured quantities by using the correspondence 
between the phase diagram and the self-bound cluster size distribution at 
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equilibrium flCampi et al.[ |2002b|) . It must be emphasized that various exper- 
iments of fragmentation of atomic nuclei^] and atomic aggregates ( |Farizon~ 
all |1998| ; |Gobet et aU , |2001| ) exhibit a percolation behaviour of the fragment 



size distribution (U-shapes, power law, and exponentials at respectively low, 
intermediate, and high energy). 



Finally, beyond the framework of simple fluids, attractive colloidal suspen- 
sions display aggregation and dynamic behaviours that could be understood 
in terms of self-bound clusters. Compared with atoms, colloidal particles have 
larger sizes, a slower diffusion, and a range of attraction much smaller than the 
hard-core diameter ( Pusey , 199 Ij ; Anderson and Lckkcrkcrker , |2002|) . These 
features contribute to form morphological structures that can be visualized 
by using microscopic imaging flSegre et aL , |2001| ). We are not aware of a mi- 



croscopic theoretical description of these clusters. What is more, percolation 
lines are frequently observed in attractive colloidal suspensions. According to 
the particular physical situation, the percolation threshold is associated with a 
sharp change in the following quantities ( |Coniglio| , |200 1| ): i) Electrical conduc- 
tivity in water- in-oil microemulsions ( |Chen et al\ |1994| ; [Weigert et at. , |1997| ). 
ii) Viscosity of triblock copolymer micellar solutions ( [Mallamace et al\ |1999| , 
[2001| ), iii) Elastic modulus of dispersions of silica spheres grafted with polymer 



chains flGrant and Russel , 1993| ), and iv) Correlation functions, measured by 
dynamic light scattering, of dispersions of silica spheres coated with stearyl 
alcohol ([Verduin and Dhont| , |1995|) . It is worth noting that in all these cases 
the percolation line passes in the vicinity of the critical point. Therefore, clus- 
tering seems to be intimately related to the dynamic properties of complex 
fluids. Clearly, more work is needed in this direction. 
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A Perfect gas of clusters model 



Let us consider a set of non-interacting clusters that do not have any volume. 
The n s indistinguishable clusters of size s form a chemical species character- 
ized by their mass ms, their chemical potential /x s and the partition function 
Qs(T,V). 

In the canonical ensemble, at constant volume and temperature, the differen- 
tial of the free energy is 

dF = ^2fx s dn s . 

s 

Furthermore, mass conservation is ensured by the condition sn s — N, or 
Y,sdn s = 0. Chemical equilibrium between clusters implies the minimization 
of the free energy. By the method of Lagrangian multipliers, we obtain 

ix s = sjx for s = l,2,...,N (A.l) 



where fx is the particle chemical potential. The fugacity z s of clusters of size s 
is then 

Zs = e ^ = e ^ = z s 

where z = z\ is the fugacity of the particles. The partition function for a given 
partition P = {ni,ri2, ■■■} of particles into non- interacting clusters is given by 

oo n 3 

Qp(T,V) = U* 1 . (A.2) 

s=l n s- 



The total partition function of a system of iV particles is then 

Q N (T, V) = J2QpKT, sn s-N) 
p 

where J2p is a sum over & h the possible partitions, whatever the number of 
particles. By moving over to the grand canonical ensemble, the grand partition 
function E(z, T, V) can be written as a function of the fugacity z: 



Z(z,T,V)=J2z N Q 



N 



N>0 

N>0 p 

oo n s 
N>0 p s=l 



By summing over N, we obtain 
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p 8=1 s=l 

p 8=1 ,b S- 

= n e (g ^ s)ns 



ln 1= ^ 



= n e 9 ^ 

8=1 

, *=i 9 < (A.3) 

The grand partition function E(z,T,V) allows us to calculate the grand po- 
tential Q = F — fiN = —kT lnH(z, T, 1/). Because f2 = —PV (see for example 
Landau and Lifshit4 |1959|) , we infer the pressure and the density in the ther- 



modynamic limit: 



1 1 00 

p P = i im ( ha) = Jim ( J2Vs(T,V)z s ) (A.4) 

z (9 In" 1 00 

p = lim ) = lim (-S^W- (A.5) 

For large enough volume, q s (T,V)/V becomes independent of V and tends 
to q s (T). In the domain of convergence of series (A.4) and (A.5), we can 
permute the limit when V — > oo and the sum over s. The pressure and density 
expression are then given by 



(3P = Y,q s (T)z s 

8=1 

oo 



8 = 1 



It is convenient to express the thermodynamic quantities as functions of the 
series n(T,z) ( [Fishery |1967a| ): 



8 = 1 

The expressions for pressure, density, compressibility, and specific heat at con- 
stant volume, become 
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PP = n(T,z) 

XT ~pdP pdzdP p 2 1 ' J 



C V = T— =-T^, = kT 2 V- 



dTv <9T 2 v> dT 2 v,„ 



Let us calculate the cluster size distribution. From now on, n s is the mean 
number of clusters of size s: 

dlnE 



dz s 

From Eq. (A.l) and Eq. (A. 3) we obtain 

n s = z s q s = q s z s . 

Define m k as the moment of order k of the cluster size distribution: 

oo 

m k = Yl skn s- 

s=l 

In this theoretical framework, thermodynamic quantities are associated with 
the moments m^. Indeed, 



1 oo 

n (T,z)M = -Y,s n qsZ s 

V s=l 



1 - n m n (T,z) 
= — > s n, = . 

We then have 



ap= m £ (^ i ) (Ag) 

,= !^ = £ (a.7) 

_ Pm 2 (T,z) _ V m 2 (T,z) 

X Vp 2 kTml(T,z) 1 j 

C^T^. (A,, 
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From Eqs. (A.6) and (A. 8) it follows thatf 37 ) 



m m 2 

PX = — = 72 

mf 



where 72 is a quantity that diverges as m 2 at the percolation threshold ( |Campi 
Tggp . 



Equations (A.6)-(A.9) clearly show that in a perfect gas of clusters model, 
geometric and thermodynamic quantities must diverge at the same points 
of the phase diagram. In particular, the second moment of the cluster size 
distribution diverges like the compressibility 



B Stability of the probabilistic Hill clusters 



To begin, we define the stability against particle emission. The separation 
energy ef s is the energy of a particle i calculated in the center of mass system of 
the cluster C s to which it belongs. A cluster is stable against particle emission 
if, for each particle i belonging to C s , we have: 



-X 



< ; • = K* + u ( r a) < 



where K* is the kinetic energy of particle i calculated in the center of mass 
system of cluster C s , that is 

= j(* - ?°') 2 , 

and V Cs = (J2iec s Vi)/ S is the velocity of the center of mass of the cluster C s . 
When the positions and velocities of the particles belonging to a given cluster 
are known, this criterion is straightforward to check. 

Let us now show that probabilistic Hill clusters are stable on average against 
particle emission in the lattice-gas model flCampi and Krivine| , |1997|) . To this 



end, we evaluate the average sign of the separation energy. For large enough 
clusters, the velocity of the center of mass can be neglected. The relative 
kinetic energy of a particle is then equal to mvf/2, where v\ is randomly given 
by the Maxwell distribution of variance l//3m. First consider the case of a 
particle interacting with only one other particle of the cluster. Its separation 
energy is then given by 

m 2 
Si = — v' 4 — e. 
2 1 



32 In a perfect gas of particles, there are only monomers: uiq = mi = = N and 
X = l/P- 
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Because the particle belongs to the cluster, we have by definition of a Hill's 
bond e.i < 0. Now, suppose the particle interacts with two particles of the 
cluster (see Fig. B.l). The separation energy of the particle is 



The probability that e« is negative is Pr{mvf/2 < 2e) = PHUi{W e )- Ori the 
other hand, the probability that the particle is linked to the cluster either by 
one or two independent bonds is given by 

2p mu (f3e) -PmuiPe). 

According to Eq. (24), we can substitute p c k for pmu- It is then very simple 
to check that 

2p ck (Pe)-p 2 ck ((3e)=p ck (2(3e). 




1) 2) 

Fig. B.l. Particle i interacts with two particles of the cluster it belongs to. 1) Particle 
i is linked to the cluster by only one Hill's bond (heavy line) with particle k, but 
interact also with particle j. 2) Particle i is linked to the cluster by two Hill's bonds 
to particles k and j . 

The generalization to more bonds between the particle and the cluster is 
straightforward. In conclusion, if a particle belongs to a probabilistic Hill 
cluster, its separation energy is negative. The condition for stability against 
particle emission is fulfilled on average. 
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